Fractal structure of a three dimensional Brownian motion on an attractive plane 
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Consider a Brownian particle in three dimensions which is attracted by a plane with a strength 
proportional to some dimensionless parameter a. We investigate the fractal spatial structure of the 
visited lattice sites in a cubic lattice by the particle around and on the attractive plane. We compute 
the fractal dimensions of the set of visited sites both in three dimensions and on the attractive plane, 
as a function of the strength of attraction a. We also investigate the scaling properties of the size 
distribution of the clusters of nearest-neighbor visited sites on the attractive plane, and compute the 
corresponding scaling exponent r as a function of a. The fractal dimension of the curves surrounding 
the clusters is also computed for different values of a, which, in the limit a — > oo, tends to that of 
the outer perimeter of planar Brownian motion i.e., the self-avoiding random walk (SAW). We find 
that all measured exponents depend significantly on the strength of attraction. 
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I. INTRODUCTION 

The laws of Brownian motion, formulated first by Ein- 
stein more than a century ago [1], have now found so 
many applications and generalizations in all quantitative 
sciences [2] . Many fractal structures in the nature can be 
derived from the sample paths of Brownian motion char- 
acterized by some appropriate fractal dimensions [3]. 

A d-dimensional Brownian motion is known to be re- 
current, i.e., the particle returns to the origin, for d <2 
and escapes to infinity for d >2. It is also known that the 
fractal (Hausdorff ) dimension of the graph of a Brownian 
motion is equal to 3/2 for d ~1, and 2 for d >2. 

The scaling limit of interfaces in various critical 2d lat- 
tice models are proven or conjectured to be described 
by the family of conformally invariant random curves 
i.e., Schramm-Loewner evolution (or SLE^) [4] which is 
driven by a Id Brownian motion of diffusivity k [5]. 

One of the most important invariance properties of pla- 
nar Brownian motion is conformal invariance. Although 
the scaling limit of 2d random walk, i.e., 2d Brownian mo- 
tion, because of self-crossing itself does not fall in the SLE 
category, variations of Brownian motion are described by 
SLE. Loop erased random walk (LERW) where loops are 
removed along the way, is one of the examples which has 
been studied by Schramm and shown that can be de- 
scribed by SLE2. 

The external perimeter of 2d random walk is also a non- 
intersecting fractal curve which can be defined by SLE. 
Verifying an earlier conjecture by Mandelbrot [3], it has 
been proven using SLE techniques [6] that the fractal di- 
mension of the Brownian perimeter is d/ =4/3, i.e, the 
same as the fractal dimension of self-avoiding random 
walk (SAW) and the external perimeter of the percola- 
tion hull. 
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In this paper, we investigate the statistical and fractal 
properties of a 3d random walker which is attracted by 
a plane. We believe that this study can provide useful 
intuitive extensions for many related physical phenomena 
including the problems with a discrete time lattice walk 
[7, 8], relaxation phenomena [9], exciton trapping [10] 
and diffusion-limited reactions [8, 11]. 



II. THE MODEL 

We consider a random walker moving along the bonds 
of a cubic lattice with the zy-plane as an attractive plane. 
The 'walker' source is considered to be the origin of the 
coordinate system. At each lattice point with z 7^ 0, 
there are six possibilities for the random walker to se- 
lect a link and move along. In our model, the random 
walker prefers walking on and near the attractive plane, 
and thus the probability that the random walker chooses 
the link which approximates it to the attractive plane is 
set to be ap, and for remaining five links is considered to 
be p, such that a > 1 (and will be called the strength of 
attraction) and p = "a^- ^"-"^ each lattice point on the 
attractive plane with z = 0, the probability that each of 
the four links on the plane to be chosen is set to be ap' 
and for two other links perpendicular to the plane is con- 
sidered to be p', where p' = ia+2 ■ ^^^ single parameter 
a in our model, controls the strength of attraction. Note 
that in the limiting case a — ?► 00 our model reduces to 
the pure 2d random walk on the plane, and for a = 1 the 
pure 3d random walk would be recovered. 
Thus there are four possible probabilities, ap' for links 
that are in the attractive plane, p' for links from the at- 
tractive plane to either of the neighboring planes, p for 
links in all of the neighboring planes or leading from them 
into the bulk, and ap for links from all the neighboring 
planes to the attractive plane. 

By detailed balance, in equilibrium at inverse temper- 
ature /3, the ratio ap/p' of the probabilities onto and 



off the attractive plane, defines an attraction energy 
/3e = ln[2a(l + 2a)/(a + 5)]. 



III. FRACTAL DIMENSION OF THE SET OF 
ALL VISITED SITES AND ITS LEVEL SET 

In the cases of random walks, systems exhibit a generic 
scale invariance, meaning that the systems can exhibit 
self-similarity and power laws without special tuning of 
parameters. This is why wc already expect that our 
model would exhibit rich fractal properties for all val- 
ues of a. 

Let us first look at the fractal spatial structure of the 3d 
attracted random walk (ARW) and its intersection with 
the attractive plane. In order to estimate the fractal di- 
mension df of the set of points visited (at least) once by 
the random walker, we examine the scaling relation be- 
tween the average number of such points Al'^^'^'^ and their 
corresponding radius of gyration Rg, i.e., M'-^'^^ ^ Rg' . 
Each ensemble averaging for M'^^'^^ (and also for M^^'^^ 
in the following) and Rg was taken over 5 x lO'* indepen- 
dent samples for a fixed number of random walk steps 
N. The measurements were done for 10'^ < N < 10^ 
with the number interval SN = 2 x 10'^. We have also 
computed the fractal dimension of the total number of 
sites on the attractive plane (i.e., M^^''^) visited by the 
random walker (in this case the corresponding radius of 
gyration is computed for all set of distinct visited sites 
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FIG. 1: (Color online) The average number of total lattice 
sites M *■ ' visited (at least) once by the attracted random 
walker (ARW) (main panel), and those M' ' on the attrac- 
tive plane (inset), as function of their average radius of gy- 
ration for two different values of the strength of attraction 
a =1.3 (■) and a =10 (a). The solid lines show the best fit 
to our data. The error bars are almost the same size as the 
symbols. 
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FIG. 2: (Color online) The fractal dimension of the set of all 
lattice points visited (at least) once by the attracted random 
walker (ARW) (■), and the set of all number of visited points 
on the attractive plane (□), as function of the strength of 
attraction a. The error bars are almost the same size as the 
symbols. 



only on the attractive plane — see Fig. 1). 
We find that the fractal dimensions have a remarkable 
continuous dependence on the parameter a. The results 
of these fractal dimensions as function of the strength of 
attraction a are illustrated in Fig. 2. As can be seen 
from figure 2, for large values of a, since the problem 
reduces to the 2d random walk on the attractive plane, 
these two fractal dimensions converge to the same value 
close to the value ~ 1.83 (this is comparable with the 
fractal dimension of the set of distinct sites visited by 
an 2d RW on a square lattice, deduced from the results 
reported in [12]). 

All error bars in this paper are estimated by using the 
standard least-squares analysis, and are almost of the 
same size as the symbols used in the figures. 
For an ideal linearly self-similar fractal of dimension df , 
one expects that the fractal dimension of the intersection 
being d'f = d/ — 1 [3] . But this is not apparently the case 
here for a ^ 1, since in our model, the attractive plane 
has disturbed the homogeneity of the probability distri- 
bution in the z-direction. Only for a ~ 1 where df = 2 



[15], we find d'. 
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IV. 



CLUSTER SIZE DISTRIBUTION ON THE 
ATTRACTIVE PLANE 



Henceforth we investigate the fractal and scaling prop- 
erties of the set of all distinct sites visited by the 3d ARW 
only on the attractive plane. Each of these sites is visited 
at least once by the 3d ARW and marked upon visiting 



(if not already). 

In this section, rather than analyzing the properties of 
the whole set, after marking all visited sites on the 
plane, we identify each cluster-site as a set of all nearest- 
neighbor visited-sites on the lattice with a specific color. 
Two typical examples of such clustering are shown in 
Fig. 3 for two different values of the strength of attrac- 
tion a = 2 and a = 10. As Fig. 3 shows, for lower 
values of a, there exist many isolated clusters of differ- 
ent scales which are accessed by the ARW only via the 
third dimension. By increasing the strength of the at- 
traction, the number of isolated clusters decreases until 
a — ^ oo for which, there will be only one large cluster on 
the attractive plane. 
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FIG. 3: (Color online) Typical samples of clusters of the vis- 
ited sites on the attractive plane by a 3d ARW of A'^ = 10^ 
shown in different colors, for q = 2 (left) and a = 10 (right). 



To examine possible scale invariance of cluster ensem- 
ble for rather small values of a, we compute the cluster 
size distribution and check whether it follows a power- 
law scaling. In the critical statistical physics, the scaling 
properties of fractal clusters can be described by the per- 
colation theory [13], where the asymptotic behavior of 
cluster distribution rig (A) near the critical point A — )■ Ac 
has the following general form 



UsiX) ^ s-^ F[{X - X,)s° 



(1) 



where a is an critical exponent, and the scaling function 
F{u) approaches to a constant value for |u| <^ 1 and 
decays rather fast for |u| ^ 1. 

We undertook simulations for several values of a to 
measure the distribution of the cluster sizes of the vis- 
ited lattice sites by the 3d ARW on the attractive plane 
(this is the probability that a visited lattice site on the 
attractive plane belongs to a cluster of size s). We gath- 
ered ensembles of a number of 5 x 10'' (for smaller a) 
and 1.5 X 10^ (for larger values of a) independent sam- 
ples of fractal patterns with marked visited-sites on the 
attractive plane. The number of the random walk steps 
was chosen to be iV = 4 x 10^ in all simulations. The 
number density Us of the clusters of size s has then been 
computed for each specific value of a by counting the 
number of clusters of size s divided by the total number 
of all clusters. 

We find that for rather small and intermediate size scale 
clusters, the distribution shows a power law behavior 
compatible with the scaling relation in Eq. (1). As can 
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FIG. 4: (Color online) Cluster size distribution exponent r 
defined in Eq. (1), as a function of the strength of attraction 
a. Inset: number density Us of clusters of the visited lattice 
sites of size s on the attractive plane for three different values 
a = 1.2, 4 and 8. The solid lines show the power-law behavior 
in the scaling region. The error bars are almost the same size 
as the symbols. 



be seen in the inset of Fig. 4, the curves for different val- 
ues of a exhibit a sharp drop-off, indicating indeed that 
they contain only small clusters. By increasing a the 
interval for scaling region decreases and a peak appears 
which signals the formation of large scale clusters. 
Our estimation of the cluster size distribution exponent 
T in the scaling region as a function of a is also shown in 
Fig. 4. One observes that the exponent r has a signifi- 
cant dependence on the strength of attraction a. 



V. FRACTAL DIMENSION OF THE CLUSTER 
BOUNDARIES ON THE ATTRACTIVE PLANE 



The remainder of this paper is dedicated to investigate 
the fractal properties of the boundaries of the visited-sites 
clusters on the attractive plane. 

Given a configuration of visited sites by the 3d ARW on 
the attractive plane, the first step is to identifying differ- 
ent clusters as outlined before. After that, the boundary 
curve of each isolated cluster has to be identified. How- 
ever the definition of interfaces and cluster boundaries on 
a square lattice can contain some ambiguities, there has 
been introduced a well-defined tie-breaking rule in [14] 
that generates non-intersecting cluster boundaries on a 
square lattice without any ambiguity. 
To define the hull for each identified cluster according to 
the algorithm defined in [14], a walker (which, of course, 
has to be distinguished from the 3d ARW) moves clock- 
wise along the edges of the dual lattice (which is also a 
square lattice) around the cluster starting from a given 
boundary edge on the cluster. The direction at each step 
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FIG. 5: (Color online) The fractal dimension of the perimeter 
of a cluster of visited sites on the attractive plane by 3d ARW, 
as a function of the strength of attraction a. Inset: the av- 
erage length of the perimeter I of a cluster versus its average 
radius of gyration r^, for two different strengths of attraction 
a = 1.2 (upper graph) and a = 16 (lower graph). The solid 
lines show the power-law behavior in the scaling region. The 
error bars are almost the same size as the symbols. 



is always chosen such that walking on the selected edge 
leaves a visited site on the right and an empty plaquctte 
on the left of the walker. If there are two possible ways 
of proceeding, the preferred direction is that to the right 
of the walker. The directions right and left are defined 
locally according to the orientation of the walker. 
According to this procedure, we have generated an en- 
semble of cluster boundary loops for several different 
strengths of attraction in the range 1.1 < a < 16. Using 

df 

the scaling relation I ^ rg , between the average length 
of the perimeter of the loops I, and their average radius 
of gyration rg , we computed the fractal dimension df of 
the cluster boundaries as a function of a. The results are 
shown in Fig. 5. 



The fractal dimension shows again a significant depen- 
dence on the strength of attraction a. In the limit a — > oo 
df converges to the value | = 1.33 which is the fractal 
dimension of the SAW i.e., the outer perimeter of the 
planar Brownian motion. 



VI. CONCLUSIONS 

In this paper, we have studied the scaling properties 
and the fractal structure of the visited lattice-sites by a 
Brownian particle in 3d which is attracted by a plane 
with the strength a. The fractal dimensions of the set of 
visited sites by the 3d random walker in both three di- 
mensions and on the attractive plane are computed which 
both converge to the same value ^ 1.83 for large a. We 
also found that size distribution of the cluster of visited 
sites by the particle on the attractive plane, has a scaling 
form characterized by an exponent that depends signifi- 
cantly on the strength of attraction. 
The fractal dimension of the surrounding loops of the 
clusters on the plane has been computed as a function 
of a. This also converges asymptotically to the expected 
value for SAW i.e., the external perimeter of a planar 
Brownian motion. 

These results need however some theoretical frame- 
work and mathematical proof. The other interesting fea- 
ture which can be investigated, is the possible confor- 
mal invariance of the cluster boundaries on the attrac- 
tive plane, which can be treated using SLE techniques 
(such study is already done only for the limiting case 
a -> oo where the problem reduces to a 2d random walk 
in the attractive plane whose boundary is described by 
SLEg /3). The fractal dimension of an SLE^ curve is given 
by df ~ 1 + k/8. In case of conformal invariance of clus- 
ter boundaries on the attractive plane in our model, they 
would be defined by a diffusivity k which depends on the 
strength of attraction. 
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